Wavevector filtering through single-layer and bilayer graphene 
with magnetic barrier structures 
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We show that the angular range of the transmission through magnetic barrier structures can be 
efficiently controlled in single-layer and bilayer graphene and this renders the structures efficient 
wavevector filters. As the number of magnetic barriers increases this range shrinks, the gaps in the 
transmission versus energy become wider, and the conductance oscillates with the Fermi energy. 
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Grapheme's electronic properties are drastically dif- 
ferent from those, say, of conventional semiconductors. 
Charge carriers in a wide single-layer graphene behave 
like " relativistic" , chiral massless particles with a "light 
speed" equal to the Fermi velocity and possess a gap- 
less, linear spectrum close to the K and K' points^. A 
perfect transmission through arbitrarily high and wide 
barriers, referred to as Klein tunneling^ is expected. In 
contrast, carriers in bilayer graphene possess a quadratic 
spectrum with zero gap. However, applying a gate allows 
the opening of a gap due to the tunnel coupling between 
the layers. This impedes the Klein tunneling and is more 
appropriate for certain applications, e. g., for improving 
the on/off ratio in carbon-based transistors. A recent 
review of the properties of graphene is given in Ref. 0. 

An angular confinement of the transmission in single- 
layer graphene was reported in Refs. andd using mag- 
netic barrier structures. However, in none of these works 
was the wavevector filtering behavior of magnetic barri- 
ers addressed. In this Letter we show how this angular 
confinement can be maximized, in single-layer and bi- 
layer graphene, and detail the resulting filtering behav- 
ior. This occurs for finite-width and <5-function barriers 
but not for complex magnetic structures in which the 
average magnetic field is zero. 

An electron in single-layer graphene, in the presence 
of a perpendicular magnetic field B(x), that varies along 
the x direction, is described by the 2x2 Hamiltonian 
Hq = vper ■ (p + eA(x)). A solution of Dirac's equation 
Hq^ = Ety, in regions where the magnetic field is not 
zero, is a linear combination of Weber functions^. For 
an electron in bilayer graphene the relevant Hamiltonian 
is a 4 x 4 matrix-! 6 - with elements H13 = H31 = t and t the 
tunnel coupling strength assumed constant t rs AOOmeV. 
Using the dimensionless units B(x) — ► B B(x), A(x) — > 
Bq£bA(x), f — ► fgf, v — > vpv, E — ► EqE, t —* Eot, with 
E = Hvf/Ib and £3 — [h/eB^] 1 / 2 , the energy spectrum 
for a homogeneous magnetic field takes the form 



For t = it gives the spectrum for single-layer graphene 6 
E = ±[2(n+ l)] 1 / 2 . The low-energy spectra for a barrier 
on single-layer and bilayer graphene are partly shown in 
the lower panel of Fig. 1 as a function of k y . For single- 
layer graphene in a region of homogeneous magnetic field 
the wave function, up to a normalization factor, is written 



in terms of the Weber functions D p (z), p 
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For bilayer graphene the corresponding wave function is 
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FIG. I: Top panel: magnetic field profile B{x) and corre- 
sponding vector potential A(x) for a series of barriers. Lower 
panel: a map in (E,k y ) space, showing the regions in which 
tunneling involving real wave vectors is allowed (in red) and 
in which is not (in yellow) for one magnetic barrier in single- 
layer and bilayer graphene. Bound states can be found in the 
yellow regions delimited by the free-particle spectra^. 
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where p ± = ( 7± - l)/2, 7± = E 12 ± (1 + i^i 2 ) 1 / 2 , and 
z = 2 x l 2 [x + k y ). In regions where the magnetic field is 
zero the solution is written as (T denotes the transpose) 

*b ± = N ± (tE, T kt ± E, kt + ik v f e ^*+*M 

(4) 

for waves propagating to the right, and as 

± = N± (tE, ±kt ± ik y ,E, -It* + ik y ) T e~ lk * x+lk y y 

(5) 

for waves propagating to the left. In this case the energy 
spectrum and wave vector kf: are 

E = Tt/2±[t 2 /4 + K ± ] 1/2 , (6) 
hi = [E 2 -k 2 y ±Etf 2 , (7) 

where K± = fc 2± + A: 2 . The factor N± = 1/iAWEk^) 1 ' 2 
can be obtained from normalization of the current. To 
evaluate the transmission through one or many barriers 
we use the transfer-matrix technique and at each inter- 
face we match the wave function and the flux 5 . Scatter- 
ing between the K and K' valleys is negligible for fields 
below 10 4 Teslal. For a single barrier we can evaluate 
the angular confinement by imposing the condition that 
the wave number k x be real for incident and transmitted 
waves. This gives 

h = [E 2 -(ky + AA 2 ] 1 / 2 , (8) 
kf = [(E±t/2) 2 -t 2 /4-(k y + A l ) 2 } 1 ^ 2 , (9) 

respectively, for a single layer and a bilayer. Further, i = 
1,2, Ai = d/2, and A 2 = -d/2. From Eqs. (8) and (9) 
and the relations k y — Ep sin 8 for a single layer and k v = 
(Ep ± Ept) 1 / 2 sin 6^ for a bilayer we can find the range 
of the angle of incidence 9 in which the transmission is 
confined. For a single barrier we obtain 

- 1 < sin 6» < 1 - d/E, (10) 

in single-layer graphene, and 

- 1 < sm6 ± < 1 -d/(E 2 ±Et) 1 / 2 . (11) 

in bilayer graphene. In single-layer graphene the angle of 
exit 0' is related to 9 by 

sin 0' = sin 8 + d/E (12) 

and in bilayer graphene by 

sin^ =sm6 ± + d/{E 2 ± Et) 1/2 . (13) 

We can also replace sin 6 by k y /k in Eqs. (10)-(13). If 
we don't use the dimensionless units, the last terms in 
Eqs. (10)-(13) are multiplied by A = i.e., 

d/E^Xd/E, d/{E 2 ±Et) 1/2 -> Xd/{E 2 ±Et) 1/2 . (14) 

From Eqs. (10) and (12) one can find the limits for 
6'. Notice that Eq. (10) implies a total reflection for 



d > 2E. We remark in passing that for the usual 2D 
Schrodinger electrons Eqs. (10) and (12) remain valid 
with E replaced by (2E) 1 / 2 . 

By controlling the ratio d/E or d/(E 2 ± Et) 1 / 2 , Eqs. 
(10) and (11) show how one can control the angular 
range of the transmission. For instance, one could 
double the width of the barrier or raise Ep, through 
the density, and/or t by doping® or gating^. However, 
one may not have a good, continuous control over these 
parameters, e.g., over the width d, which may be limited 
by fabrication techniques. It is therefore worth searching 
for other ways to control these parameters. With this in 
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FIG. 2: Contour plot of the transmission through one mag- 
netic barrier of width d — 4£b in (a) and through four barriers 
of width d — £b in (b) for single- layer graphene. (c), (d): Re- 
spectively as in (a) and (d) for bilayer graphene with t = 15. 



mind we investigate the effect of increasing the number 
of barriers on the transmission. The region in which 
tunneling, involving real wavevectors k x , is allowed is 
shaded in red in Fig. 1. For n barriers the angular range 
is given by Eqs. (10) and (11) upon changing d to nd. 
From Eqs. (10) and (11) we can evaluate the minimum 
energy required for tunneling. For n barriers of width d 
and 8 = we have E = nd/2 for single-layer graphene 
and E = T{t/2) ± (t/2)[l + n 2 d 2 /t 2 \ x l 2 for bilayer 
graphene. For t nd we can rewrite this minimum as 
E w ±n 2 d 2 /At. This shows clearly that we can control 
the tunneling by increasing the number of barriers. 
Typically the height of the magnetic barriers is in the 
range 0.1-lTesla. For B — IT the magnetic length is 
lb = 25.6 nm and the energy scale Eq = 25.7 meV for a 
typical vp = 10 6 m/s. In Fig. 2 we show a contour plot 
of the transmission, in (E, k y ) space, for one barrier of 
width d = 4\£b and for four magnetic barriers of width 
d = Ib- As seen using four barriers instead of one results 
in sharp resonances shown in the lower part of the right 
panels for both single-layer and bilayer graphene. This 
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FIG. 3: (a) Contour plot of the transmission in [E, 8) space 
for four barriers of width d = Ib with L — 10£b in bilayer 
graphene. The inset shows the result for one barrier with 
d — Ib- (b) As in (a) but for single-layer graphene. 



holds for any multibarrier structure n > 2. As Eq. (10) 



experiments than the transmission. For very low temper- 
atures, such that f(E) - f(E + eV) ps -eVS(E - E F ), 
with f(E) the Fermi-Dirac function and eV the 
source-to-drain voltage drop, the standard formula 
for G is G = G Q J^ 2 /2 T(E F ,E F sme)cos9de, with 

Go = 2e 2 EpL y /irh and L y the length of the struc- 
ture along the y direction. In Fig. 5(a) we plot the 
conductance of a four-barrier structure in single-layer 
graphene and in Fig. 5(b) in bilayer graphene with 
d = £b and two values of the interbarrier distance, 
L = 1Q£b and L — 20£b, and t = 15. The insets show 
the corresponding results for a single barrier of width 
d = 4£b- As shown, in line with the transmission results 
of Fig. 3 those for the conductance show oscillations in 
the multibarrier structure that are absent in the single 
barrier. In summary, we showed that the angular range 
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FIG. 4: Transmission probability for different number of 
magnetic barriers in bilayer graphene (a) and single-layer 
graphene (b) with the parameters of Fig. 3. 



shows, the angular range of the transmission shrinks by 
increasing the number of barriers. We show that in Fig. 
3(a) for bilayer graphene and in Fig. 3(b) for single-layer 
graphene. Had we used one barrier of width d = A£b, the 
shrinking of the shown angular range would be the same 
but the resonaces at the lower left end of the plot would 
be absent. This can be seen more clearly in Fig. 4 where 
we plot the transmission versus the energy for 8 = — n/4 
and different number of barriers. Notice that the energy 
gap widens with n. Having seen the transmission 
resonances, one may wonder to what extent they survive 
the averaging over the Fermi surface that determines 
the conductance G which is more easily accessed in 



FIG. 5: Conductance through four magnetic barriers, of 
width d — Ib, in single-layer graphene (a) and in bilayer 
graphene (b) for two values of L and t — 15. The upper in- 
sets show the result through one magnetic barrier of width 
d = Mb. 



of the transmission through magnetic barrier structures 
can be efficiently controlled by increasing the number 
of barriers in single-layer and bilayer graphene. This 
renders the structures efficient wavevector filters, leads 
to transmission resonances and conductance oscillations, 
and widens the gaps in the transmission as a function of 
the energy. 
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